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Abstract. We reveal cosmic star-formation history obscured by dust using deep infrared observa- 
tion with the AKARI. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 
/""N | 15, 18, and 24/xm) by the AKARI satellite allows us to estimate restframe 8/im and 12/im lumi- 

nosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in 
previous work. We found that restframe 8jUm (0.38 < z < 2.2), 12^im (0.15 < z < 1.16), and total 
infrared (TIR) luminosity functions (LFs) (0.2 < z < 1.6) constructed from the AKARI NEP deep 
<*S !■ data, show a continuous and strong evolution toward higher redshift. In terms of cosmic infrared 

luminosity density (£2/#), which was obtained by integrating analytic fits to the LFs, we found a 
good agreement with previous work at z < 1.2, with Q.ir°^ ( 1 + z) AA± 1 °. When we separate contri- 
butions to Q,jr by LIRGs and ULIRGs, we found more IR luminous sources are increasingly more 
' important at higher redshift. We found that the ULIRG (LIRG) contribution increases by a factor of 

10 (1.8) from z=0.35 to z=1.4. 
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Introduction. Revealing the cosmic star formation history is one of the major goals 
of the observational astronomy. However, UV/optical estimation only provides us with 
a lower limit of the star formation rate (SFR) due to the obscuration by dust. A straight- 
forward way to overcome this problem is to observe in infrared, which can capture the 
star formation activity invisible in the UV. The superb sensitivities of recently launched 
Spitzer and AKARI satellites can revolutionize the field. 

However, most of the Spitzer work relied on a large extrapolation from 24/im flux to 
estimate the 8, 12/im or total infrared (TIR) luminosity, due to the limited number of 
mid-IR filters. AKARI has continuous filter coverage across the mid-IR wavelengths, 
thus, allows us to estimate mid-IR luminosity without using a large ^-correction based 
on the SED models, eliminating the largest uncertainty in previous work. By taking 
advantage of this, we present the restframe 8, 12/im TIR LFs, and thereby the cosmic 
star formation history derived from these using the AKARI NEP-Deep data. 

Data & Analysis. The AKARI has observed the NEP deep field (0.4 deg 2 ) in 9 filters 
(N2,N3,N4,S7,S9W,SU,L15,Lim and L24) to the depths of 14.2, 11.0, 8.0, 48, 58, 
71, 117, 121 and 275juJy (5a)dl. This region is also observed in BVRi'z' (Subaru), 
u' (CFHT), FUV,NUV (GALEX), and J,Ks (KPN02m), with which we computed 
photo-z with =0.043. Objects which are better fit with a QSO template are removed 
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FIGURE 1. (left) Restframe 8/im LFs. The blue diamonds, purple triangles, red squares, and orange 
crosses show the 8/im LFs at 0.38 < z < 0.58,0.65 < z < 0.90,1.1 < z < 1.4, and 1.8 < z < 2.2, 
respectively. The dotted lines show analytical fits with a double-power law. Vertical arrows show the 
8/im luminosity corresponding to the flux limit at the central redshift in each redshift bin. Overplotted 
are Babbedge et al. %M in the pink dash-dotted lines, Caputi et al. 121 in the cyan dash-dotted lines, 
and Huang et al. [6] in the green dash-dotted lines. AGNs are excluded from the sample, (middle) 
Restframe 12/im LFs. The blue diamonds, purple triangles, and red squares show the 12/im LFs at 
0.15 < z < 0.35,0.38 < z < 0.62, and 0.84 < z < 1.16, respectively. Overplotted are Perez-Gonzalez 
et al. UH at z=0.3,0.5 and 0.9 in the cyan dash-dotted lines, and Rush, Malkan, & Spinoglio 1 12] at z=0 
in the green dash-dotted lines, (right) TIR LFs. 



from the analysis. We compute LFs using the 1/V ma x method. Data are used to 5 a with 
completeness correction. Errors of the LFs are from 1000 realization of Monte Carlo 
simulation. 



8/J.m LF. Monochromatic 8/im luminosity (Ls^ m ) is known to correlate well with 
the TIR luminosity OIlaL especially for star- forming galaxies because the rest-frame 
8/im flux are dominated by prominent PAH features such as at 6.2, 7.7 and 8.6 /im. The 
left panel of FigJT] shows a strong evoltuion of 8/im LFs. Overplotted previous work had 
to rely on SED models to estimate L^ m from the Spitzer S^/im in the MIR wavelengths 
where SED modeling is difficult due to the complicated PAH emissions. Here, AKARI's 
mid-IR bands are advantageous in directly observing redshifted restframe 8/im flux in 
one of the AKARI's filters, leading to more reliable measurement of 8/im LFs without 
uncertainty from the SED modeling. 



12jim LF. 12/im luminosity (L\2^m) represents mid-IR continuum, and known to 



correlate closely with TIR luminosity II 1 ill - The middle panel of FigJT] shows a strong 



evoltuion of 12/im LFs. Here the agreement with previous work is better because (i) 
12/im continuum is easier to be modeled, and (ii) the Spitzer also captures restframe 
12/im in SiA^im at z=l. 

TIR LF. Lastly, we show the TIR LFs in the right panel of Fig[TJ We used Lagache, 
Dole, & Puget JSC's SED templates to fit the photometry using the AKARI bands 
at >6/im (57,S9W,S11,L15,L18W and L24). The TIR LFs show a strong evolution 
compared to local LFs. At 0.25 < z < 1.3, L* TIR evolves as (1 +z) 4 - 1±0 - 4 . 




FIGURE 2. Evolution of TIR luminosity density based on TIR LFs (red circles), 8/xm LFs (stars), and 
12/xm LFs (filled triangles). The blue open squares and orange filled squares are for LIRG and ULIRGs 
only, also based on our Ljir LFs. Overplotted dot-dashed lines are estimates from the literature: Le Floc'h 
et al. H, Magnelli et al. [10|] , Perez-Gonzalez et al. |0, Caputi et al. J2], and Babbedge et al. J2] are 
in cyan, yellow, green, navy, and pink, respectively. The purple dash-dotted line shows UV estimate by 
Schirninovich et al. [13]. The pink dashed line shows the total estimate of IR (TIR LF) and UV lfl3ll . 



Cosmic star formation history. We fit LFs in Fig{T]with a double-power law, then 
integrate to estimate total infrared luminosit y d ensity at various z. The restframe 8 
and 12/im LFs are converted to Ljir using [ 11, 0] before integration. The resulting 
evolution of the TIR density is shown in Figf2l The right axis shows the star formation 
density assuming Kennicutt |7|]. We obtain Q.ir(z) (1 +z) 4 ' 4±L0 . Comparison to £luy 
fl 1 3M suggests that Qtir explains 70% of 0. to tai at z=0.25, and that by z=1.3, 90% of 
the cosmic SFD is explained by the infrared. This implies that Q.tir provides good 
approximation of the Q. tota i at z > 1 . 

In Figl2l we also show the contributions to Qtir from LIRGs and ULIRGs. From 
z=0.35 to z=1.4, Q./R by LIRGs increases by a factor of ~1.6, and Q. 1R by ULIRGs 
increases by a factor of ~10. More details are in Goto et al. [H]. 



Spatially-Resolved Spectroscopy of an E+A (post-starburst) System. We per- 
formed a spatially-resolved medium resolution long-slit spectroscopy of a nearby E+A 
(post-starburst) galaxy system with FOCAS/Subaru 0]. This E+A galaxy has an obvi- 
ous companion galaxy 14kpc in front (FiglH left) with the velocity difference of 61.8 
km/s. 

We found that H5 equivalent width (EW) of the E+A galaxy is greater than 7A galaxy 
wide (8.5 kpc) with no significant spatial variation. We detected a rotational velocity in 
the companion galaxy of >175km/s. The progenitor of the companion may have been 
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FIGURE 3. (left) The SDSS g, r, /-composite image of the J1613+5103. The long-slit positions are 
overlayed. The E+A galaxy is to the right (west), with bluer colour. The companion galaxy is to the left 
(east), (right) H5 EW is plotted against D4000. The diamonds and triangles are for the E+A core/north 
spectra, respectively. The squares and crosses are for the companion galaxy's core/north spectra. Gray 
lines are population synthesis models with 5-100% delta burst population added to the lOG-year-old 
exponentially-decaying (T=lGyr) underlying stellar population. SalpeterlMF and metallicity of Z = 0.008 
are assumed. On the models, burst ages of 0.1, 0.25, 0.5 and 2 Gyr are marked with the filled circles. 



a rotationally-supported, but yet passive SO galaxy. The age of the E+A galaxy after 
quenching the star formation is estimated to be 100-500Myr, with its centre having 
slightly younger stellar population. The companion galaxy is estimated to have older 
stellar population of >2 Gyrs of age with no significant spatial variation (FiglH right). 

These findings are inconsistent with a simple picture where the dynamical interaction 
creates infall of the gas reservoir that causes the central starburst/post-starburst. Instead, 
our results present an important example where the galaxy-galaxy interaction can trigger 
a galaxy-wide post-starburst phenomena. 
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